UCI-TR-2003-6 



Superweakly Interacting Massive Particles 
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We investigate a new class of dark matter: superweakly-interacting massive particles (super- 
WIMPs). As with conventional WIMPs, super WIMPs appear in well-motivated particle theories 
with naturally the correct relic density. In contrast to WIMPs, however, super WIMPs are impossible 
to detect in all conventional dark matter searches. We consider the concrete examples of gravitino 
and graviton cold dark matter in models with supersymmetry and universal extra dimensions, re- 
spectively, and show that superWIMP dark matter satisfies stringent constraints from Big Bang 
nucleosynthesis and the cosmic microwave background. 
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There is ample evidence that luminous matter makes 
up only a small fraction of all matter in the universe. Re- 
sults from the Wilkinson Microwave Anisotropy Probe, 
combined with other data, constrain the non-baryonic 
dark matter density to £1dm = 0.23 ± 0.04 |lj, far in ex- 
cess of the luminous matter density. We therefore live in 
interesting times: while the amount of dark matter is be- 
coming precisely known, its identity remains a mystery. 

WIMPs, weakly- interacting massive particles with 
weak-scale masses, are particularly attractive dark mat- 
ter candidates. WIMPs have several virtues. First, 
their appearance in particle physics theories is indepen- 
dently motivated by the problem of electroweak symme- 
try breaking. Second, given standard cosmological as- 
sumptions, their thermal relic abundance is naturally 
that required for dark matter. Third, the requirement 
that WIMPs annihilate efficiently enough to give the de- 
sired relic density generically implies that WIMP-matter 
interactions are strong enough for dark matter to be dis- 
covered in current or near future experiments. 

Here we consider a new class of non-baryonic cold dark 
matter: superweakly-interacting massive particles (su- 
perWIMPs or S WIMPs). As with WIMPs, superWIMPs 
appear in well-motivated theoretical frameworks, such as 
supersymmetry and extra dimensions, and their (non- 
thermal) relic density is also naturally in the desired 
range. In contrast to conventional WIMPs, however, they 
interact superweakly and so evade all direct and indirect 
dark matter detection experiments proposed to date. 

For concreteness, we consider two specific super- 
WIMPs: gravitinos in supersymmetric theories, and 
Kaluza-Klein (KK) gravitons in theories with extra di- 
mensions. Gravitino and graviton superWIMPs share 
many features, and we investigate them in parallel. 

For gravitino superWIMPs, we consider supergravity, 
where the gravitino G and all standard model (SM) su- 
perpartners have weak-scale masses. Assuming i?-parity 
conservation, the lightest supersymmetric particle (LSP) 
is stable. In supergravity, the LSP is usually assumed 
to be a SM superpartner. Neutralino LSPs are excel- 
lent WIMP candidates, giving the desired thermal relic 



density for masses of 50 GeV to 2 TeV, depending on 
Higgsino content. In contrast, here we assume a G LSP. 
The gravitinos considered here couple gravitationally and 
form cold dark matter, in contrast to the case in low-scale 
supersymmetry breaking models where light gravitinos 
couple more strongly and form warm dark matter. 

We consider also the possibility of graviton dark mat- 
ter in universal extra dimensions (UED), in which gravity 
and all SM fields propagate 0] • We focus on D = 5 space- 
time dimensions with coordinates x AI = {x^,y). The 
fifth dimension is compactified on the orbifold S 1 /Z2, 
where S 1 is a circle of radius R, and Z2 corresponds 
to y — > —y. Unwanted massless fields are removed 
by requiring suitable transformations under y — > —y. 
For example, the 5D gauge field Vm{%, y) transforms 
as Vp(x,y) -> Vn(x,y) and V 5 (x,y) -* -V 5 (x,y) under 
y — > —y, which preserves V®(x) and removes V®(x). Sim- 
ilar choices remove half of the fermionic degrees of free- 
dom, producing chiral 4D fcrmions, and preserve the 4D 
graviton /i°„(x) while removing h^ 5 (x) and h^^x). The 
gravi-scalar h® 5 (x) remains; we assume that some other 
physics stabilizes this mode and generates a mass for it. 

The orbifold compactification breaks KK number con- 
servation, but preserves KK-parity. KK particles must 
therefore be produced in pairs, and current bounds re- 
quire only R^ 1 > 200 GeV 0,0] ■ KK-parity conserva- 
tion also makes the lightest KK particle (LKP) stable 
and a dark matter candidate. For ~ TeV, weakly- 
interacting KK particles have thermal relic densities con- 
sistent with observations 0. In particular, B 1 , the first 
KK partner of the U(l) gauge boson, has been shown to 
be a viable WIMP dark matter candidate, with promis- 
ing prospects for direct detection 0, and also indi- 
rect detection in anti-matter searches |5j, neutrino tele- 
scopes 000], and gamma ray detectors . 

As in the case of supersymmetry, however, the lightest 
partner need not be a SM partner. In UED, the LKP 
could be G , the first KK partner of the graviton. G 1 is, 
in fact, perhaps the most natural LKP candidate — ra- 
diative contributions to KK masses, typically positive 0, 
are negligible for G 1 . G 1 couplings are also gravitational, 
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and so highly suppressed. 

Gravitinos and gravitons therefore naturally emerge as 
super WIMP candidates: stable massive particles with su- 
perweak interactions. Their weak gravitational interac- 
tions imply that they play no role in the thermal history 
of the early universe. (We assume inflation followed by 
reheating to a temperature low enough to avoid regener- 
ating large numbers of super WIMPs.) Thus, if the next 
lightest supersymmetric particle (NLSP) or next lightest 
KK particle (NLKP) is weakly-interacting, it freezes out 
with a relic density of the desired magnitude. Much later, 
however, these WIMPs then decay to superWIMPs; as 
the WIMP and superWIMP masses are similar, the su- 
perWIMP then inherits the desired relic density. 

Unlike WIMPs, however, superWIMPs are impossi- 
ble to discover directly, and their annihilation rate is 
so suppressed that they also escape all indirect detec- 
tion experiments. At the same time, unlike superheavy 
dark matter candidates with only gravitational interac- 
tions [I3> superWIMPs inherit the desired relic density 
from a thermal abundance and arise from accessible elec- 
troweak physics. At colliders, WIMP decays to super- 
WIMPs will occur long after the WIMP leaves the detec- 
tor. If the NLSP or NLKP is neutral, its metastability 
will have no observable consequences. The discovery of a 
seemingly stable but charged NLSP or NLKP may, how- 
ever, provide a strong hint for superWIMP dark matter. 

We now investigate constraints on and alternative sig- 
nals of superWIMP dark matter scenarios. The observ- 
able consequences of superWIMPs must rely on the de- 
cays of WIMPs to superWIMPs on cosmological time 
scales [llj. The NLSP or NLKP may be any SM partner. 
In supergravity, the lightest SM superpartner is often the 
Bino £?, the superpartner of the hypercharge gauge bo- 
son. In the minimal UED scenario |2Ll|, the lightest SM 
KK mode is often B 1 . Motivated by these results, we 
now consider specific scenarios in which decays to super- 
WIMPs are typically accompanied by photons, and we 
consider the impact of electromagnetic cascades. 

In the supersymmetric photon superWIMP sce- 
nario, NLSP decay is governed by the coupling 
-drfi^^ P \^BF vpi where F is the U(l) field 
strength, and M* = {^Gn)^ 1 ' 2 ~ 2.4 x 10 18 GeV is 
the reduced Planck scale. The NLSP decay width is 
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In models with low-scale supersymmetry breaking, 
and the gravitino couples dominantly 



ttlq -C m B , 

through its ±h spin components. In the high-scale su- 
persymmetry breaking scenarios considered here, how- 
ever, the couplings of the ±| spin polarizations are of 
the same order and must be kept in deriving Eq. (|T)l. 

The properties of gravitons in UED scenarios may 
be determined straightforwardly; details will be pre- 
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FIG. 1: Lifetimes for B — > G-y (left) and B 1 -> G lr y (right) 
for Am = mwiMP — itiswimp and mswiMP = 0.1 TeV (long 
dashed), 0.3 TeV (short dashed), and 1 TeV (solid). 



sented elsewhere [l2j. Graviton superWIMPs couple to 
B 1 through &Gl„(-F°> J -PF lv p + \r)^ F^F^), where 
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F™„ = dpEZ - d u B*. The G 1 B 1 B° vertex is identical to 
the G°B°B° vertex, and the longitudinal component of 
the massive B 1 plays no role. The NLKP decay width is 
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The B and B 1 lifetimes are given in Fig. ^ In both 
cases, in the limit Am = m WIMP - toswimp "C toswimp, 
the WIMP lifetime is proportional to (Am) -3 and is in- 
dependent of the overall mass scale. 

As we will see, the most relevant bounds constrain the 
total energy released in photons in WIMP decay, or more 
precisely, e 7 Ky, where e 7 is the energy of the photons 
when created and Y y = n 7 /n^ G is the number density 
of photons from WIMP decay normalized to the num- 
ber density of background photons n^ G = 2£(3)T 3 /-7r 2 , 
where T is the temperature during WIMP decay. In the 
superWIMP scenario, WIMPs decay essentially at rest, 
and so e 1 = (m^ IMP - m swiMp)/ ( 2m wiMp)- Since a su- 
perWIMP is produced in association with each photon, 
Y" 7 = Yswimp, and the photon abundance is given by 



Iswimp si 3.0 x 1CT 12 



TeV 
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^SWIMP 



0.23 
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Predicted values for £ 7 1swimp are shown in Fig. [21 

As evident in Fig. ^ WIMP decays occur long after 
Big Bang nucleosynthesis (BBN), and so may, in prin- 
ciple, destroy the successful BBN predictions for light 
element abundances. The energy of photons produced in 
late decays is rapidly redistributed through scattering off 
background photons, 77rq — > e + e~, and inverse Comp- 
ton scattering [l3l Il4| . As a result, the constraints of 
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FIG. 2: The photon energy release £ 7 Yswimp for various 
mswiMP in TeV in the gravitino (left) and graviton (right) 
superWIMP scenarios. We fix Oswimp = 0.23; £ 7 Yswimp 
scales linearly with fiswiMP- BBN constraints exclude the 
shaded regions [113 ]: consistency of the CMB with a black- 
body spectrum excludes regions above the CMB contours. 



BBN are, to an excellent approximation, independent of 
the initial energy distribution of injected photons, and 
constrain only the total energy release. 

Detailed analysis [ljj, demanding consistent predic- 
tions for deuterium, 3 He, 4 He, 6 Li, and 7 Li, excludes 
the region of parameter space shown in Fig. [3 The BBN 
constraint is weak for early decays: at early times, the 
universe is hot and the initial photon spectrum is rapidly 
softened, leaving few high energy photons above thresh- 
old to modify the light element abundances. We find 
that BBN excludes some of the relevant parameter re- 
gion, but not all. In particular, for relatively short-lived 
WIMPs with r < 10 7 s and weak-scale masses, the re- 
quirement that superWIMPs form all of the dark matter 
is consistent with the constraints from BBN. 

The cosmic microwave background (CMB) also im- 
poses constraints |l5l Il6j . The injection of energy in the 
form of photons may distort the CMB from the observed 
black-body spectrum. Before redshifts of z ~ 10 7 , elastic 
Compton scattering, bremsstrahlung eX — > eXj (with 
X an ion), and double Compton scattering e~j — > e-77 
effectively thermalize injected energy. After z ~ 10 7 , 
however, the photon number-changing interactions be- 
come ineffective, and the photon spectrum relaxes only 
to a Bose-Einstein distribution with chemical potential fi. 
After z ~ 10 5 , even Compton scattering becomes ineffec- 
tive, and deviations from the black-body spectrum may 
be parametrized by the Sunyaev-Zeldovich y parameter. 

As with the BBN constraints, bounds from the CMB 
are largely independent of the injected energy spectrum, 
depending primarily on the total energy release. The 
bounds on £ 7 Y"swimp scale linearly with the bounds on 
/i and y. We update the analysis of Ref. 0] to include 
the latest results \fj,\ < 9 x 10" 5 and|w| < 1.2 x 1CT 5 [l7| . 
with baryon density fi^/i 2 ~ 0.022 |l|, where h ~ 0.71 is 



the normalized Hubble expansion rate. These bounds ex- 
clude energy releases above the CMB contours in Fig. [2] 
Remarkably, the CMB constraints are now so precise that 
they supersede BBN constraints for decay times r ~ 10 7 s 
and t > 10 10 s. Nevertheless, regions of superWIMP pa- 
rameter space, including regions with weak-scale masses 
and mass splittings, remain viable. 

Finally, for highly degenerate WIMP-SWIMP pairs, 
WIMPs decay very late to soft photons. The photon 
spectrum is not thermalizcd and may produce observ- 
able peaks in the diffuse photon spectrum. The present 
differential flux of photons from WIMP decay is 



dt iV W iMp(i) 
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where to ~ 13.7 Gyr is the age of the universe Q], 

N W iMp(t) = -^wiMP e_ *^ rvvlMP » where ^wimp is the num- 
ber of WIMPs at freeze-out, and Vq is the present volume 
of the universe. The diffuse photon flux is a sensitive 
probe only when WIMPs decay in the matter-dominated 
era. We may then take 1 + z = (to/t) 2 ' s , and 
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where F(a) = a 1/2 e a , a = (E 7 /e 7 )(io/ r wiMp) 2/3 - 
F(a) is maximal at a = 3~ 2 / 3 , and so, for twimp < to, 
the differential photon flux reaches its maximal value at 
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= 680 keV 



GeV 
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for gravitinos, and an energy 1.4 times smaller for gravi- 
tons. The Am dependence follows from the redshifting 



of photons created with energy Am by 1 + z oc r, 
(Am) 2 . For high degeneracies, 
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and so the maximal flux is 

— (£") = 1.5 cm" 2 s- 1 sr- 1 MeV" 
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for gravitinos, and a factor of 1.4 larger for gravitons. 

Representative photon energy spectra are shown in 
Fig. |3J Also shown are the measured diffuse fluxes from 
the observatories HE AO, OSSE, and COMPTEL [l|, 
determined from observed fluxes by subtracting known 
point sources. The observed spectra fall rapidly with en- 
ergy and so severely constrain the relatively hard photon 
spectra predicted by Am < 10 GeV. Note that we have 
not included photon interactions which soften the photon 
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FIG. 3: Diffuse photon fluxes (solid) for mswiMP = 1 TeV, 
^swimp = 0.23, and Am = 1 GeV (solid) and 10 GeV (long 
dashed) , and upper bounds from observations (short dashed) . 




Diffuse photon 

I k I I 777 ^ 

0.1 0.5 1.0 5.010.0 50.00. 1 
m smM p/(f!swiMp/0.23) (TeV) 



Diffuse photon 

I I 77 

0.5 1.0 5.010.0 50.0 

swimp/(^swimp/0-^3) (TeV) 



FIG. 4: Regions of the (mswiMP, Am) plane excluded by 
BBN, CMB, and diffuse photon constraints. The shaded re- 
gions and the regions below the CMB contours are excluded. 



spectrum for Am > few GeV such effects can only 
enlarge the allowed parameter space discussed below. 

In Fig. we compile the constraints discussed above 
and show the allowed regions of the (toswimp, Am) plane 
for the photon gravitino and graviton superWIMP dark 
matter scenarios. The BBN and CMB constraints are 
as discussed above. An additional region is excluded by 
the requirement that the diffuse photon flux never exceed 
the observed flux by 2a for any energy. Although these 
data exclude some of the parameter space, the most well- 
motivated region with mswiMP, Am ~ 100 GeV to 1 TeV 
remains an outstanding possibility. 

In conclusion, we find that superWIMPs provide a 
qualitatively novel possibility for particle dark matter. 
Such particles appear in the form of gravitinos and gravi- 
tons in theories with supersymmetry and extra dimen- 
sions, and they naturally inherit the desired relic density 
from late-decaying weakly-interacting NLSPs or NLKPs. 

SuperWIMPs satisfy existing constraints from BBN 



and CMB and evade all conventional dark matter ex- 
periments. On the surface, we have apparently com- 
mitted Pauli's "ultimate sin" by proposing a solution to 
the dark matter problem that has no observable conse- 
quences. However, improvements in the measurements 
discussed above may uncover anomalies. For example, a 
detailed study of the 15 keV to 10 MeV diffuse photon 
flux by INTEGRAL is currently underway 0. A pro- 
nounced bump in this spectrum could provide a striking 
signal of superWIMP dark matter. Finally, we note that 
neutrino, charged lepton, weak gauge boson, and Higgs 
boson NLSPs and NLKPs are all viable from the point of 
view of preserving the naturalness of the desired relic den- 
sity. Some of these scenarios may be severely constrained 
by bounds on hadronic showers, but the remaining sce- 
narios will have qualitatively different, and possibly very 
interesting, observational consequences. 

We thank J. Cohn, T. Han, T. Moroi, and M. White 
for enlightening conversations. 
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